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I.  INTRODUCTION 


Lee.  Hobbs  and  .Atkinson  [1]  have  developed  a  computer  program  TRUCK  3.1  to  find  the 
<.1\  aainic  response  of  a  vehicle  subjected  to  different  types  of  loads  including  such  intensive 
loads  a.'  blast  waves.  The  code  is  intended  to  yield  gross  motions  of  the  vehicle  body,  and 
of  tire^.  a.xles.  shelters,  and  racks  relative  to  the  vehicle  body.  Large  motions,  including  the 
o\erturning  of  the  vehicle,  are  permitted.  However,  each  individual  element  like  an  axle, 
sheltei  or  rack  is  treated  as  a  rigid  body. 

The  purposes  of  this  study  are:  (i)  to  check  the  validity  of  the  governing  equations  of 
nioiiuii  Used  in  the  TRUCK  code  that  model  the  rigid  body  motions  of  the  vehicle  assembly 
and  its  components,  (ii)  to  find  out  if  the  equations  have  been  correctly  coded  and  (iii)  to 
investigate  whether  or  not  the  TRUCK  code  models  effectively  the  frictional  force  between 
the  tires  and  the  ground. 

The  third  objective  is  achieved  by  analyzing  a  few  sample  problems  and  comparing  their 
computed  solution  with  either  analytical  solutions  or  those  obtained  by  using  the  computer 
code  .-\DI.\.-\. 


U.  VEHICLE  MODEL 

The  vehicle  and  its  primary  components  whose\esponse  can  be  studied  by  using  the 
TRUCK  code, are  shown  in  Fig.  l._  Each  component  is  characterized  by  its  ma^s,  center  of 
gravity  position,  inertias  and  products  of  inertias.  'Not  explicitly  shown  in  Fig.  1  are  springs 
connecting  the  vehicle  body  to  the  axles.''One  might  think  that  the  mechanical  properties  of 
these  springs  will  affect  the  motion  of  axles,  racks  etc.  relative  to  the  vehicle  body  but  will  not 
affect  the  overall  motions  of  the  assembly.  But  in  the  TRUCK/"  code,  the  primary  coordinate 
system  is  a  body  axis  system  attached  to  the  vehicle  body  with  its  origin  at  the  cente’-  of 
gravity  of  the  entire  system.  Therefore,  the  mechanical  properties  of  springs  connecting  the 
vehicle  body  to  the  axles  influence  the  overall  motions  of  the  vehicle.  However,  the  motion 
of  the  overall  center  of  gravity  of  a  frictionless  system  is  not  affected  by  the  mechanical 
properties  of  these  springs.  The  right  handed  rotations  <^,  9,  and  rp  about  the  i,  y,  and 
;-axes,  shown  in  Fig.  1,  are  referred  to  as  pitch,  yaw,  and  roll,  respectively.  Translational 
motions  in  the  x,  y,  and  z  directions  will  be  referred  to  as  sideslip,  heave,  and  fore-and-aft 
translation,  respectively.  If  desired,  the  motions  of  the  axles,  shelters,  and  racks  relative  to 
the  vehicle  body  can  also  be  studied. 

Each  axle  is  assigned  two  degrees  of  freedom,  namely  its  roll  and  heave  relative  to  the 
vehicle  body.  A  shelter  is  considered  to  have  four  degrees  of  freedom  -  its  roll,  heave,  pitch 
and  fore-and-aft  translation  relative  to  the  vehicle  body.  A  rack  has  six  degrees  of  freedom  - 
its  roll,  sideslip,  heave,  pitch,  yaw  and  fore-and-aft  translation  relative  to  the  vehicle  body. 

1  he  1  RUCK  code  has  the  option  to  consider  the  entire  system  as  a  rigid  body  with 
only  three  degrees  of  freedom  involving  vehicle  roll,  heave  and  sideslip  only.  However,  this 
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problem  still  requires  data  on  the  tires  necessary  to  define  the  tire/ground  interaction  forces. 

III.  FRICTIONAL  MODELS  IN  TRUCK  AND  ADINA  CODES 

In  the  TRUCK  code,  the  tire-ground  interaction  is  presumed  to  be  based  upon  Cou¬ 
lomb's  law  of  friction.  The  normal  force  on  the  tire  depends  upon  the  normal  deflection  and 
normal  velocity  of  the  tire.  For  this  purpose  the  tire  is  modeled  as  a  non-linear  spring  and  a 
non-linear  dashpot  connected  in  series.  The  magnitude  of  the  tangential  force  depends  upon 
whether  the  tire  is  sliding  or  not.  When  there  is  no  sliding,  the  tangential  force  depends 
upon  the  tire  tangential  deflection  and  velocity  and  the  spring  and  damping  characteristics 
of  the  tire.  The  tire  starts  to  slide  when  the  tangential  spring  force  equals  the  tire  normal 
force  times  the  coefficient  of  friction.  The  tangential  spring  force  during  sliding  equals  the 
tire  normal  force  multiplied  by  the  coefficient  of  friction  and  is  in  a  direction  opposite  to  the 
tangential  deflection. 

We  note  that  in  the  TRUCK  code  the  static  and  kinetic  coefficients  of  friction  are  given 
the  same  value. 

The  tangential  velocities  and  therefore  tangential  displacements  of  front  tires  arc  treated 
differently  from  those  of  other  tires.  This  is  because  of  the  assumption  that  the  front  tires 
are  free  to  roll  in  their  planes.  Therefore,  for  a  front  tire,  the  component  of  the  tangential 
velocity  which  is  normal  to  the  intersection  of  the  wheel  plane  and  the  ground  plane  is  used 
rather  than  the  total  tangential  velocity. 

The  possibility  of  a  tire  leaving  the  ground  and  coming  in  contact  again  is  allowed. 
For  a  tire  that  is  off  the  ground,  tire-ground  interaction  forces  are  set  equal  to  zero  and  its 
tangential  deflection  is  also  set  equal  to  zero. 

In  the  AE84  version  of  the  ADINA  [2]  code  one  of  the  two  contacting  bodies  must  be 
deformable.  For  two-dimensional  problems,  the  code  sets  the  frictional  force  between  the 
two  bodies  equal  to  that  required  to  equilibrate  the  entire  system  provided  that  it  does  not 
exceed  the  limiting  value  of  the  coefficient  of  friction  times  the  normal  force  between  them. 
In  the  latter  case,  the  frictional  force  is  set  equal  to  the  limiting  value  and  sliding  is  permitted 
between  the  two  bodies. 

In  ADINA  code  too,  there  is  no  distinction  made  between  the  static  and  kinetic  coeffi¬ 
cients  of  friction. 

Herein  we  will  study  those  problems  in  w’hich  the  vehicle  undergoes  planar  motions. 

FVU  COMPUTATION  AND  DISCUSSION  OF  RESULTS 

.After  having  verified  that  equations  governing  the  motion  of  the  vehicle  assembly  and 
its  components  are  modeled  correctly  in  the  TRUCK  code  we  proceeded  to  find  out  if  the 
code  predicts  w^ell  the  response  of  an  idealized  system  subjected  to  simple  loading.  The 


\ehK'l(‘  coii^iiiereci  hoiireforth  is  sliown  in  Fig.  2  wherein  are  also  listed  values  of  \arious 
paranieteis.  In  order  to  mininuze  the  axle  roll  relati'e  to  the  vehicle  body  the  axles  were 
siven  arbitrarily  higli  ■.  .ues  of  L-.  When  computing  the  numerical  solution  by  using  the 
I  Kl  t'K  rxide.  tlie  ..routine  that  gives  aerodynamic  loading  on  the  vehicle  was  bypassed. 

.\'  a  lir''  te^i  problem,  we  aiijdiixl  a  force  F  (as  a  function  of  time  t)  to  the  rear  of 
the  vehicle  ^o  a^  to  ])ro(!uce  fore-and-aft  translation  only,  and  set  the  coefficient  of  friction 
//  l.iet  ween  the  tire-,  and  the  ground  equal  to  zero.  The  force  is 

f  =  4.45  t  k  MN  =  k  lb 

where  k  i-  a  unit  vector  in  the  c-direetioii.  Thus  the  problem  reduces  to  finding  the  motion 
of  a  rigid  block  of  mass  rn  (having  units  of  kg  (Ib  s^/inj)  sliding  on  a  smooth  plane  and 
subjected  to  the  external  force  /  .  Its  acceleration  a,  velocity  r  and  displacement  u  relative 
to  the  earth-tixed  axes  are  gi\'en  by 
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The  comptited  and  analytical  .solutions  for  the  fore-and-aft  translation  are  sliown  in  Fig.  3. 
The  two  solutions  agreed  with  each  other.  .Also  the  computed  normal  tire  forces  matched 
exactly  with  those  found  analytically.  This  exercise  ensured  that  governing  equations  are 
integrated  correctly  at  least  for  this  simple  loading. 

Subsequently  wr  cl’.aiiged  the  value  of  /r  from  0  to  0.8.  The  maximum  frictional  force 
between  tfie  vehicle  and  tfie  ground  for  the  assumed  weight  of  the  vehicle  assembly  is  151.1  kN 
(33.970  Ibl.  Thus,  theoretically,  the  vehicle  should  not  slide  till  the  applied  force  exceeds 
151.1  k.\  (33.970  lb).  However,  the  frictional  force  between  the  tires  and  the  ground  need 
not  equal  the  applied  force  in  magnitude  but  it  will  act  in  the  opposite  direction.  Fig.  4 
shows  tlie  magnitude  of  the  computed  frictional  force  and  the  applied  force.  As  dictated  by 
Coulomb's  friction  law.  the  vehicle  did  not  slide.  Ihe  oscillations  of  the  computed  frictional 
force  about  the  applied  force  are  due  to  the  modeling  of  the  system  as  springs,  dashpots. 
and  connecting  rigid  masses. 

In  order  to  see  how  the  assumed  tangential  tire  stiffness  will  affect  the  computed  solu¬ 
tion.  we  repeated  the  above  calculation  with  two  other  values  of  the  tangential  tire  stiffness. 
The  results  depicted  in  Fig.  5  reveal  that  the  computed  frictional  force  depends  upon  the 
tangential  tire  stiffness. 

WV'  sulisequcnt ly  investigated  the  effect  of  the  size  of  the  time  increment  upon  the 
conqinted  solution.  For  tirtie  increment  equal  to  one-tenth  of  its  previous  value,  the  computed 
frictional  force  tor  three  values  of  Ay  is  shown  in  Fig.  b.  3  he  results  shown  in  Figs.  5 
and  6  indicate  that  the  maffriitiide  of  the  lime-increment  to  be  used  in  the  TRTCK  code 
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•  Damping  constants  for  tires  =  0. 

•  Tires  are  modeled  as  linear  springs  with  spring  constants  of  17.5  GN/m  (10*  Ib/in)  for 
both  normal  and  tangential  displacements. 

•  Time  increment  =  lOOps. 

•  Damping  constants  for  springs  between  axle  and  body  =  0. 

•  Spring  constants  for  springs  between  axle  and  body  =  1.75  GN/m  (10^  Ib/in). 

•  Linear  springs,  used  to  model  the  tire  stiffness  normal  to  the  ground,  can  slide  on  the 

ground. 


Figure  2.  Vehicle  Model  Studied. 


ELAPSED  TIME  (10'“  s) 


Figure  3.  Com  ,r  ;on  of  the  Analytical  Solution  with  that  Obtained  by  Using  the  TRUCK 
Code  (No  Fr;ct  ■  j  Force  between  Tires  and  Ground). 


FORCE  (lb) 


elapsed  time  (10-''  s) 


Figure  4.  Comparison  of  the  Applied  Force  and  the  Frictional  Force  Computed  with  the 
TRUCK  Code. 
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FORCE  (kN) 


FORC  E  (lb) 


FORCE  (kN) 


FORCE  (lb) 


FORCE  (kN) 


ELAPSED  TIME  (lO--*  s) 

Figure  7.  Comparison  of  Frictional  Force  Computed  with  the  TRUCK  Code  and  Applied 
Loads  when  Tires  Slide. 

depends  strongly  upon  the  tangential  stiffness  of  the  tire.  The  rather  odd  result  for  k^  — 
1.75  TN/m  =  10^°  Ib/in  could  be  due  to  improper  choice  of  the  time  increment  At.  We 
have  not  investigated  this  any  further.  We  note  that  our  values  for  fc,  are  hypothetical  and 
need  not  correspond  to  a  real  tire. 

In  Fig.  7  are  plotted  the  values  of  the  computed  frictional  force  and  the  applied  force 
when  the  coefficient  of  friction  was  lowered  from  0.8  to  0.05.  The  limiting  frictional  force 
equals  9.448  kN  (2,124  lb).  Because  of  the  symmetries  in  the  assumed  vehicle  model,  each 
wheel  should  provide  2.362  kN  (531  lb)  of  frictional  resisting  force  in  the  limiting  case. 
However,  because  of  the  special  treatment  rendered  to  front  tires,  the  TRUCK  code  gave 
the  limiting  frictional  force  to  be  4.724  kN  (1,062  lb)  when  the  rear  wheels  started  sliding. 
This  is  exactly  half  of  the  theoretical  value  2md  the  difference  between  the  two  is  due  to  the 
way  frictional  forces  are  computed  for  the  front  and  the  rear  tires.  Note  that  the  front  tires 
did  not  slide  and  apparently  contributed  nothing  to  the  frictional  force. 

When  solving  the  same  problem  with  the  ADINA  code,  the  tires  were  assumed  to  have 
negligible  radius.  Thus,  the  entire  bottom  surface  was  in  contact  with  the  flat  ground.  The 
ADINA  code  gave  values  of  frictional  force  which  were  in  agreement  with  Coulomb’s  friction 
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Figure  8.  Discrete  Models  Used  in  the  TRUCK  and  ADINA  Codes. 

law.  That  is,  when  the  externally  applied  force  was  less  than  the  limiting  value  of  the 
frictional  force  between  the  vehicle  and  the  ground,  the  code  gave  frictional  force  equal  in 
magnitude  and  opposite  in  direction  to  the  applied  force.  When  the  applied  force  exceeded 
the  magnitude  of  the  limiting  frictional  force,  the  computed  frictional  force  equalled  the 
limiting  frictional  force. 

A  possible  explanation  for  the  difference  in  the  computed  frictional  force  with  the  two 
codes  is  that  the  equivalent  discrete  models,  shown  in  Fig.  8,  used  in  the  two  codes  are  quite 
distinct.  In  the  ADINA  code  the  rigid  block  was  modeled  as  a  deformable  body  with  very 
high  value  of  Young’s  modulus.  An  attempt  to  model  tires  by  equivalent  truss  elements  in 
the  ADINA  code  did  not  give  satisfactory  results.  As  already  discussed  modelling  of  tires 
as  very  strong  springs  in  the  TRUCK  code  necessitates  the  use  of  an  extremely  small  time 
increment. 

Note  that  in  the  discrete  model  used  with  the  TRUCK  code,  there  are  springs  connecting 
the  axles  to  the  body  of  the  vehicle.  In  the  ADINA  code,  the  entire  vehicle  was  modeled 
as  a  rigid  body.  The  modelling  of  the  vehicle  by  springs  and  dashpots  with  interconnecting 
masses  contributes  to  the  oscillations  of  the  frictional  force  computed  with  the  TRUCK  code. 
An  attempt  was  made  to  investigate  the  effect  of  the  stiffness  of  the  springs  between  the 
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axles  and  the  vehicle  body  upon  the  computed  frictional  force.  The  foregoing  results  were 
obtained  by  taking  =  1.75  GN/m  =  10"  ib/in  and  At  =  lOO^s.  When  we  increased  the 
value  of  k^  to  1.75  T.\/m  (10^*^  Ib/in).  numerical  instabilities  were  observed  for  At  =  lOOfis. 
These  instabilities  subsided  for  a  while  with  At  =  10/us  and  At  =  1/iS  but  appeared  again  at 
a  later  instant.  It  seems  that  one  will  need  smaller  time-step  size  to  pursue  this  case  further. 
V\e  should  add  that  during  the  time  interval  for  which  stable  results  were  obtaiined  with  the 
TRUCK  code,  the  computed  frictional  force  essentially  matched  with  the  applied  force  and 
also  with  that  computed  with  the  ADINA  code. 

Even  though  there  were  no  dampers  included  in  either  of  the  discrete  models,  the 
amplitude  of  the  oscillations  of  the  frictional  force  computed  wdth  the  TRUCK  code  seemed 
to  die  out  possibly  due  to  the  dissipation  introduced  by  the  numerical  method  used  to 
integrate  the  equations  of  motion. 

In  the  loading  situation  studied  above  with  the  TRUCK  code,  the  rear  tires  slid  whereas 
the  front  ones  did  not.  One  situation  where  all  four  tires  are  likely  to  slide  simultaneously 
is  when  the  force  is  applied  to  a  side  of  the  truck.  We  tried  to  simulate  this  by  applying  a 
force 

F  =  4.45ti  MN  =  lO^tilb 

at  the  center  of  gravity  of  the  truck  assembly.  Values  of  all  other  parameters  are  given  in 
Fig.  2.  The  computed  frictional  forces  for  /i  =  0.8  and  /i  =  0.05  are  plotted  in  Fig.  9.  As 
expected,  all  tires  slid  simultaneously  in  the  latter  case  and  the  computed  frictional  force 
did  eventually  stabilize  near  the  limiting  value  of  9.448  kN  (2,124  lb). 

To  see  how  well  the  TRUCK  code  predicts  the  rotational  response  of  a  vehicle,  we 
changed  the  loading  in  the  vehicle  model  of  Fig.  2  to  zero  applied  force  but  applied  a  yawing 
moment  M  given  by 

M  =  1.13  tj  MN-m  =  10^  tj  lb-in 

Results  were  computed  for  F,  =  17.5  GN/m  (10®  Ib/in),  /u  =  0.8;  kf  =  17.5  GN/m  (10®  Ib/in), 
/I  =  0.1;  and  F,  =  175  GN/m  (10®  Ib/in),  /i  =  0.8.  For  /i  =  0.8  and  0.1,  the  limiting 
frictional  moments  equal  345.4  kN-m  (30.57x10®  lb-in)  and  43.16  kN-m  (3.82x10®  lb-in), 
respectively.  The  computed  values  of  the  moment  caused  by  frictional  forces  are  plotted 
in  Fig.  10.  Results  for  —  17.5  GN/m  (10®  Ib/in)  and  /i  —  0.8  essentially  agreed  with 
those  computed  for  /i  =  0.1.  However,  when  fc*  was  changed  from  17.5  GN/m  (10®  Ib/in)  to 
175  GN/m  (10®  Ib/in)  with  /i  kept  fixed  at  0.8,  the  computed  values  of  the  frictional  moment 
varied  widely.  These  are  listed  in  Table  1. 

These  simple  exercises  indicate  that  with  a  suitable  choice  of  the  values  of  the  tangential 
tire  stiffness  and  the  time  step  size,  the  TRUCK  code  gives  good  results.  It  will  be  desirable 
to  build  into  the  code  an  algorithm  that  optimizes  the  time  increment  and  checks  for  the 
stability  of  the  computed  results. 

We  should  add  that  neither  the  applied  force  nor  the  applied  moment  produced  over¬ 
turning  of  the  vehicle  in  the  simple  problems  studied. 

For  3-dimensional  problems,  the  AE84  version  of  the  ADINA  code  models  contact 
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FORCE  (lb) 


tt  =  10'“’  sec..  TRUCK  CODE  USED 


ELAPSED  TIME  (10-“  s) 


Figure  9.  Comparison  of  Computed  Frictional  Force  and  Applied  Force  when  the  External 
Force  Is  Applied  to  the  Side  of  a  Truck. 
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FORCE  (kN) 


MOMENT  AEOUT  Y  AXIS  (K)''  lb-in 


Figure  10.  Comparison  of  the  Applied  Moment  and  Frictional  Moment  Computed  with  the 
TRUCK  Code. 
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MOMENT  ABOUT  Y-AXIS  (kN- 


conditions  as  either  frictionless  or  sticking  i.e.  with  infinite  friction.  Once  the  additional 
capability  to  account  for  finite  non-zero  frictional  forces  at  the  contact  surfaces  has  been 
included  in  the  ADLXA  code,  it  should  be  possible  to  solve  problems  involving  overturning 
of  a  vehicle  with  both  the  TRUCK  and  ADIIS’A  codes  and  compare  the  results.  We  note 
that  all  of  the  problems  studied  above  were  2-dimensional  and  the  ADIN.A.  code  modeled 
frictional  effects  properly. 


V.  CONCLUSIONS 

The  equations  of  motion  used  in  the  TRUCK  code  were  found  to  be  correct  and  accu¬ 
rately  coded  in  the  computer  program.  The  code  gave  acceptable  values  of  frictional  forces 
between  the  tires  and  the  ground. 

For  the  problems  studied,  the  ADINA  code  gave  frictional  forces  which  agreed  with 
those  obtained  by  using  Coulomb's  law. 

Problems  involving  overturning  of  a  vehicle  were  not  solved  with  either  the  TRUCK  or 
the  ADIX.A  code. 

The  results  indicate  that  the  magnitude  of  the  time-increment  to  be  used  in  the  TRUCK 
code  depends  strongly  upon  the  tangential  stiffness  of  the  tire. 
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